Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes

We have developed the formalism necessary to employ the discontinuous-Galerkin approach in general-relativistic hydrodynamics. The formalism is firstly presented in a general 4-dimensional setting and then specialized to the case of spherical symmetry within a 3+1 splitting of spacetime. As a direct application, we have constructed a one-dimensional code, EDGES, which has been used to asses the viability of these methods via a series of tests involving highly relativistic flows in strong gravity. Our results show that discontinuous Galerkin methods are able not only to handle strong relativistic shock waves but, at the same time, to attain very high orders of accuracy and exponential convergence rates in smooth regions of the flow. Given these promising prospects and their affinity with a pseudospectral solution of the Einstein equations, discontinuous Galerkin methods could represent a new paradigm for the accurate numerical modelling in relativistic astrophysics.Comment: 24 pages, 19 figures. Small changes; matches version to appear in PR

Multimessenger Parameter Estimation of GW170817

We combine gravitational wave (GW) and electromagnetic (EM) data to perform a Bayesian parameter estimation of the binary neutron star (NS) merger GW170817. The EM likelihood is constructed from a fit to a large number of numerical relativity simulations which we combine with a lower bound on the mass of the remnant's accretion disk inferred from the modeling of the EM light curve. In comparison with previous works, our analysis yields a more precise determination of the tidal deformability of the binary, for which the EM data provide a lower bound, and of the mass ratio of the binary, with the EM data favoring a smaller mass asymmetry. The 90\% credible interval for the areal radius of a $1.4\ M_\odot$ NS is found to be $12.2^{+1.0}_{-0.8} \pm 0.2\ {\rm km}$ (statistical and systematic uncertainties).Comment: 7 pages, 3 figures, accepted to the EPJA Topical Issue: The first Neutron Star Merger Observation - Implications for Nuclear Physic

Universality and intermittency in relativistic turbulent flows of a hot plasma

With the aim of determining the statistical properties of relativistic turbulence and unveiling novel and non-classical features, we resent the results of direct numerical simulations of driven turbulence in an ultrarelativistic hot plasma using high-order numerical schemes. We study the statistical properties of flows with average Mach number ranging from $\sim 0.4$ to $\sim 1.7$ and with average Lorentz factors up to $\sim 1.7$. We find that flow quantities, such as the energy density or the local Lorentz factor, show large spatial variance even in the subsonic case as compressibility is enhanced by relativistic effects. The velocity field is highly intermittent, but its power-spectrum is found to be in good agreement with the predictions of the classical theory of Kolmogorov. Overall, our results indicate that relativistic effects are able to significantly enhance the intermittency of the flow and affect the high-order statistics of the velocity field, while leaving unchanged the low-order statistics, which instead appear to be universal and in good agreement with the classical Kolmogorov theory. To the best of our knowledge, these are the most accurate simulations of driven relativistic turbulence to date.Comment: 5 pages, 4 figures. Minor changes to match the version accepted on ApJ

Universality and intermittency in relativistic turbulent flows of a hot gas

With the aim of determining the statistical properties of relativistic turbulence and unveiling novel and non-classical features, we present the results of direct numerical simulations of driven turbulence in an ultrarelativistic hot plasma using high-order numerical schemes. We study the statistical properties of flows with average Mach number ranging from $\sim 0.4$ to $\sim 1.7$ and with average Lorentz factors up to $\sim 1.7$. We find that flow quantities, such as the energy density or the local Lorentz factor, show large spatial variance even in the subsonic case as compressibility is enhanced by relativistic effects. The velocity field is highly intermittent, but its power-spectrum is found to be in good agreement with the predictions of the classical theory of Kolmogorov.Comment: Talk given at the ASTRONUM2012 conference on the 25th of June 201

High-Order Numerical-Relativity Simulations of Binary Neutron Stars

We report simulations of the inspiral and merger of binary neutron stars performed with \texttt{WhiskyTHC}, the first of a new generation of numerical relativity codes employing higher than second-order methods for both the spacetime and the hydrodynamic evolution. We find that the use of higher-order schemes improves substantially the quality of the gravitational waveforms extracted from the simulations when compared to those computed using traditional second-order schemes. The reduced de-phasing and the faster convergence rate allow us to estimate the phase evolution of the gravitational waves emitted, as well as the magnitude of finite-resolution effects, without the need of phase- or time-alignments or rescalings of the waves, as sometimes done in other works. Furthermore, by using an additional unpublished simulation at very high resolution, we confirm the robustness of our high convergence order of $3.2$.Comment: Submitted for the ASTRONUM-2014 proceedings. Includes a previously unpublished high-resolution simulatio